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1. 


INTRODUCTION 


This  report  describes  the  application  of  a  new  type 
of  eigenfunction  expansion  to  the  analysis  of  sound 
propagation  in  a  two-dimensional  waveguide.  In  the 
present  discussion  one  boundary  of  the  waveguide  is  flat, 
the  other  boundary  is  uneven,  and  the  waveguide  is  filled 
with  inhomogeneous,  lossless  fluic*.  The  expansions  are 
designed  to  accommodate  the  conditions  that  the  pressure 
vanish  at  the  flat  boundary  and  that  the  normal  component 
of  fluid  velocity  vanish  at  the  uneven  boundary. 

The  method  of  expansion  is  novel  in  that  two  essen¬ 
tially  independent  functions  are  expanded  simultaneously. 
Expansions  of  this  kind  provide  a  complete  representation 
of  the  exact  pressure  and  velocity  fields  at  the  bound¬ 
aries  of  the  waveguide  as  well  as  in  the  fluid.  Conditions 
at  the  uneven  boundary  can  be  satisfied  because  the  eigen¬ 
value  problem  that  generates  the  expansion  functions  has  a 
boundary  condition  that  contains  the  eigenvalue.  This 
boundary  condition  is  the  key  to  our  analysis.  It  causes 
the  eigenvalues  to  be  complex  even  though  the  waveguide 
is  lossless. 

The  coefficients  in  the  field  expansions  vary  along 
the  direction  of  propagation.  These  coefficients  satisfy 
a  linear  system  of  first  order  ordinary  differential 
equations  that  are  almost  identical  to  coupled  mode 
equations  for  waveguides  that  have  parallel  plane  bound¬ 
aries.  The  coefficients  are  also  discontinuous  at 
positions  associated  with  line  sources  in  the  waveguide. 

To  emphasize  the  treatment  of  the  uneven  boundary, 
we  assume  in  the  present  discussion  that  the  fluid 
density  is  a  continuous  function  of  depth  alone.  This 
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restriction  is  not  fundamental  to  our  methods,  however. 

In  our  next  report  we  shall  treat  lateral  variations  in 
density  as  well. 

The  present  report  is  organized  as  follows.  First, 
we  discuss  the  general  eigenvalue  problem  and  state  the 
principal  eigenfunction  expansion  theorem  associated 
with  it  (Section  2) .  Although  the  eigenfunctions  are 
not  orthogonal  in  any  classical  sense,  there  is  a  simple 
expression  in  closed  form  for  the  expansion  coefficients. 
Next,  we  begin  discussion  of  the  waveguide  problem  with  a 
description  of  its  mathematical  formulation  (Section  3) . 
Eigenvalue  problems  of  the  kind  reviewed  at  the  outset  of 
this  report  originate  in  an  asymptotic  analysis  of  self- 
sustaining  oscillations  in  the  waveguide  (Section  4) . 

These  problems  generate  the  eigenfunctions  that  we  use 
to  represent  the  exact  sound  fields  in  the  waveguide. 
Finally,  we  derive  a  linear  system  of  first  order 
ordinary  differential  equations  along  with  jump  conditions 
that  govern  the  coefficients  in  the  eigenfunction 
expansions  of  the  exact  sound  fields  (Section  5) . 

2.  SURVEY  OF  THE  EIGENVALUE  PROBLEM 

The  eigenvalue  problem  associated  with  self-sustaining 
oscillations  in  the  waveguide  is  a  two-point  boundary  value 
problem  for  a  second  order  ordinary  differential  equation 
on  a  finite  interval.  The  eigenvalue  appears  quadratically 
in  the  differential  equation  and  linearly  in  one  boundary 
condition.  The  differential  equation  has  the  form 


P  (z)  — 

*  y 

—  —  <Mz) 

+ 

' 

K2( z)  -  k2 

dz 

,  p(z)  dz 

2 


on  the  open  interval  -H  <  z  <  0.  We  require  that  the 
functions  p(z)  and  K{z)  take  only  real  values,  and  that 
p(z)  >0.  In  addition,  we  assume  that  dp(z)/dz  and  K(z) 
are  continuous  on  the  closed  interval  -H^z^O.  The 
boundary  conditions  are 


—  <f>  ( — H )  +  Hik  <M-H)  =  0  , 

(2) 

dz 

<j>  (0)  =  o  . 

(3) 

The  parameter  H  may  take  any  real  value. 

A  detailed  study  of  the  spectrum  of  eigenvalues  k 
will  be  presented  in  the  next  report.  Here  we  briefly 
note  some  of  the  main  features  of  the  spectrum.  First 
of  all,  it  is  discrete.  We  may  list  the  eigenvalues  as 
kn,  where  the  index  n  takes  any  value  among  the  positive 
and  negative  whole  numbers  ±1,  ±2,  and  so  on.  Second, 
all  but  at  most  a  finite  number  of  the  eigenvalues  are 
purely  imaginary.  The  distribution  of  those  eigenvalues 
having  nonvanishing  real  parts  is  symmetric  with  respect 
to  the  imaginary  axis  in  the  complex  k  plane.  These 
eigenvalues  are  bounded  in  magnitude  by  the  maximum 
value  that  the  function  |K(z)|  attains  on  the  closed 
interval. 

|k  |  <  max  | K ( z ) |  if  Re (k  )  ^  0  .  (4) 

n  -H<z<° 

If  H  M  then  the  imaginary  parts  of  these  eigenvalues 
do  not  vanish.  This  is  clear  from  the  identity 
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H  |  <T(-H)  | 


Im(k  )  = 
n 


,2  .  .  P(-H) 

<j>  (z)  |  _ dz 


if  Re  (kn)  /  0  ,  (5) 


-H 


P  (z) 


and  the  fact  that  4>n(~H)  /  0;  for  condition  (2)  would 
force  the  eigenfunction  to  vanish  identically  if  4>n  ( — H ) 
were  zero.  Finally,  the  purely  imaginary  eigenvalues 
kn  tend  to  be  spaced  almost  periodically  at  intervals 
of  7t/H  along  the  imaginary  axis  in  the  complex  k  plane 
when  the  magnitude  of  kn  is  large. 

The  eigenvalues  can  be  listed  in  a  way  that  makes 
description  of  some  of  these  features  easy.  For  example, 
the  eigenvalues  with  nonvanishing  real  parts  can  be 
grouped  together.  There  is  an  even  number,  say  2v,  of 
such  eigenvalues,  and  the  listing  can  be  constructed  so 
that  Re(kn>  ft  o  it  and  only  if  |n|  <_  v.  In  addition,  it 
is  convenient  to  order  the  eigenvalues  of  this  kind  so 
that 


k  =  -k*  if  Re(k  )  ft  0  , 
-n  n  n 


(6) 


where  the  asterisk  denotes  the  complex  conjugate.  The 
purely  imaginary  eigenvalues  can  be  listed  with  regard 
to  the  trend  among  them  eventually  to  form  two  infinite, 
almost-periodic  sequences: 


k  =  i[0  +  (n-O.5)n]/H  +  0(  1/n)  as  n  -*■  +°°  , 
n 


kn  =  i[8+(n  +  0.5)iT]/H  +  0(  1/n)  as  n  . 


(7) 

(8) 
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The  angle  6  is  determined  by  the  parameter  H. 

0  =  arctan(H)  where  jej  <  rr/2  .  (9) 

A  simple  analysis  explains  the  asymptotic  behavior 
of  the  purely  imaginary  eigenvalues  and  of  the  corresponding 
eigenfunctions.  For  large  |k  |  we  introduce  approximate 

a  a  n 

solutions  k  and  $  (z)  such  that 
n  Yn 


(10) 

(11) 


The  lowest  order  approximation  to  problem  (1-3)  in  terms 

/V  A 

of  kn  and  4>n(z)  has  the  same  form  as  if  p(z)  were  constant 
and  K ( z)  vanished  identically: 


JL_  l  (z)  -  k  2  $  (z)  =  0  ,  (12) 

di7  n  n  n 


lin(-H)  +  iiiMn(«H)  =  0  ,  (13) 


$  (0)  =  0  .  (14) 

n 


Thus, 


A 

k 

n 


i  < 


(15) 


4>n  ( z )  =  Ansin(<nz)  ,  ( 16) 
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where  the  amplitude  An  /  0  is  arbitrary,  and 

Kn=[0+(n-O.5)7T]/H  if  n  =  1,  2,  ...  ,  (17 

<n=[0+(n  +  O.5)TT]/H  if  n  =  -1,  -2,  -  .  (18 

Equations  (10-11)  and  (15-18)  therefore  account  for  the 
approximate  periodicity  of  the  purely  imaginary  eigen¬ 
values,  and  they  also  indicate  that  the  corresponding 
eigenfunctions  are  basically  sinusoidal.  These  features 
suggest  that  there  might  be  enough  redundancy  among  the 
eigenfunctions  for  them  to  be  used  to  represent  two 
functions  at  the  same  time.  For  example,  if  H  =  0  then 
the  spectrum  can  be  divided  into  two  disjoint  parts,  and 
the  eigenfunctions  associated  with  each  part  can  be 
used  to  construct  a  unique  representation  of  a  single 
function.  In  general  such  a  division  appears  not  to  be 
convenient.  We  shall  see,  however,  that  the  collection 
of  eigenfunctions  can  always  be  easily  used  as  a  whole 
to  construct  unique  representations  of  two  functions 
simultaneously. 

Orthogonality  usually  plays  a  role  in  theories  of 
eigenfunction  expansions.  There  is  no  orthogonality 
relation  in  any  classical  sense  associated  with  the 
eigenfunctions  of  problem  (1-3)  if  H  /  0,  however.  In 
fact,  any  one  of  these  eigenfunctions  can  be  expressed 
infinitely  many  ways  in  terras  of  the  others.  Neverthe¬ 
less,  in  the  theory  of  these  eigenfunctions  there  is  a 
relation  involving  both  the  eigenfunctions  and  also  the 
eigenvalues  which  is  somewhat  similar  to  orthogonality. 

In  conjunction  with  certain  convergence  properties  of 
the  eigenfunction  expansions,  this  relation  can  be  used 
to  calculate  the  expansion  coefficients.  It  may  there¬ 
fore  be  regarded  as  a  natural  substitute  for  the  lacking 


orthogonality  principle.  We  state  it  in  the  next  theorem. 


Theorem. 

Let  k  and  k  be  eigenvalues  of  problem  (1-3) ,  and 
m  n 

let  4>  (z)  and  4>n(z)  be  the  corresponding  eigenfunctions. 

(a)  If  k  ^  k  then 
m  n 


rO 

(k  +  k  )  <J>  (z)<|>(z)£-^-dz  -  Hi  <P  (-H)<t> 

mn  mnp(z)  m  n 


)  =  0  .  (19) 


(b)  For  all  but  at  most  a  finite  number  of  indices  n, 


2  k 

n 


f° 

4>2(z) 

n 

-H 


P  ( ~H) 
P  (z) 


dz 


Hi*  (~H) 
n 


f  0  . 


(20) 


To  derive  relation  (19)  let  us  first  express  eq.(l) 

in  operator  form  L(<j>»k)  =  0.  We  then  obtain  the  following 

standard  identity  by  partial  integration  of  the  difference 

(<(>  /p)LU  ,k  )  -  (<b  /p)L(<)>  ,k  ): 
m  n  n  n  mm 


rO 

(k2  -  k2)  <}>  (z)  $  (z)  _i _ dz 

m  n  Lh  m  n  P  (z) 


<Pn(z) 


— - —  <t>  ( z )  -  4>  (z)  — - —  <P  (z) 

p(z)  dz  m  m  p  ( z )  dz  n 


z  =  0 


z  =  -H 


(21) 


Application  of  boundary  conditions  (2-3)  simplifies  the 
right  side  of  this  identity. 
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T 


i 

t 

\ 

■i 

\ 


I 


>m(z)*n(z) 


dz 


-H 


P  (z) 


=  (k  -  k  )  H  i  $  (-H)  <|>  (-H) 
m  n  m  n 


P  (-H) 


(22) 


Then  division  by  the  factor  (k^  -  kn)  /  0  yields  eq.(19). 

To  derive  relation  (20)  let  us  first  observe  from 
the  asymptotic  approximations  (10-11)  that 


2  k 


-H 


f^(z)  — — ^LL  dz  -Hi  <'2(-H) 
P  (  2)  n 


=  2  k 


-H 


$2(z)  dz  -  Hi4>2(-H)  +  O  MA2|/n 
n  in 


(23) 


as  n  -*■  +<*.  Evaluation  of  the  right  side  of  this  approx¬ 
imation  by  use  of  solutions  (15-18)  is  straightforward. 


rO 


2  k 


2  2, 


n 


'( z )  d  z  - 


H  i  <j>  (-H) 
n 


k  HA 
n  n 


(24) 


-H 


Equations  (23-24)  imply  eq.(20)  since  |  k  |  -*•  00  as  n  -*•  ±°°. 

We  conclude  this  section  with  a  statement  of  the 
principal  eigenfunction  expansion  theorem  associated  with 
problem  (1-3).  For  a  pair  of  essentially  independent 
functions  f(z)  and  g(z)  on  the  closed  interval  -H<_z  <_0, 
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this  theorem  asserts  the  existence  and  the  uniqueness  of 
eigenfunction  expansions  of  the  form 


f(z)  =  l  an4>n(z)  ,  (25) 

n  =  -oo 

oo 

g  (z)  =  l  anikn<J>n(z)  •  (26> 

n  =  -oo 

(The  accent  indicates  that  a  sum  has  no  n  =  0  term, 
in  accord  with  our  way  of  listing  the  solutions  of  the 
eigenvalue  problem.)  The  coefficients  an  appearing  in 
these  expansions  are  determined  jointly  by  the  functions 
f(z)  and  g(z).  The  expansions  are  valid  simultaneously; 
the  same  coefficients  an  appear  in  both  of  them.  We  may 
therefore  regard  the  eigenfunction  expansions  (25-26)  as 
a  representation  of  the  complex  vector  valued  function 
(f(z),g(z))  in  terms  of  the  nonorthogonal  basis  elements 
(<)>n(z)  ,  i  kn  <t>n(z) )  . 

In  the  statement  of  this  theorem  we  shall  assume 
that  relation  (20)  is  valid  for  all  n.  This  assumption 
is  not  fundamental  to  our  theory.  Indeed,  in  the  next 
report  we  shall  discuss  some  cases  that  violate  it,  and 
we  shall  state  and  prove  the  eigenfunction  expansion 
theorems  that  hold  for  these  cases.  To  simplify  the 
introductory  discussion  of  our  expansion  methods,  however, 
we  have  chosen  to  avoid  these  complications  here. 

The  Principal  Eigenfunction  Expansion  Theorem. 

Suppose  that  the  function  f(z)  has  a  continuous  third 
derivative,  and  that  the  function  g(z)  has  a  continuous 
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second  derivative.  Su  pose  that  these  functions  also 
satisfy  the  boundary  conditions 


d 

—  f  (-H)  +  Hg(-H)  =  0  , 

(27) 

dz 

f(0)  =  0  , 

(28) 

g(0)  =  o  . 

(29) 

If  relation  (20)  is  valid  for  all  n  in  the  realization 
of  problem  (1-3)  that  generates  the  eigenfunctions,  then 
there  exist  unique  eigenfunction  expansions  of  the  form 
(25-26)  for  the  functions  f(z)  and  g(z).  Furthermore, 
the  expansion  for  f(z)  can  be  differentiated  twice,  and 
the  expansion  for  g(z)  can  be  differentiated  once. 

Parts  (a)  and  (b)  of  the  theorem  make  these  statements 
precise . 

(a)  Existence  of  the  eigenfunction  expansions. 

Define  the  coefficients  an  by  means  of  the  formula 


a 

n 


f  (z)  -  i  g  (z)]<J>  (z) 
n 


P(-H) 
P  (z) 


dz  -  H  i  f  (-H)  <}>  ( -H) 
n 


/0 


2  k 


n 


-H 


9  P(-H)  .  2 

<{/(z) - dz  -  H  i  <|>  ( -H) 

n  P (z)  n 


for  n  =  ±1,  ±2,  and  so  on  . 


(30) 


(The  assumption  about  relation  (20)  prevents  division  by 
zero  in  this  formula.) 
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Define  two  sequences  of  functions  {f  (z)}  and  {g  (z)} 

N  N 

in  terms  of  these  coefficients  by  means  of  the  symmetric 
partial  sums 

N 

fN(z)  =  l  an*n(2)  ' 

n  =  -N 

N 

VZ)  =  J  anikn*n(2)  • 

n  =  -N 

Then  the  following  sequences  of  functions  converge 
uniformly  to  the  indicated  limits  on  the  closed  interval 
-H_<z_<0  as  N  -*•  oo ; 

j  j 

f  (z)  -*■  f  (z)  ,  g  (z)  •>  g(z)  and  —  f  (z)  ■+  — f  (z)  .  (33) 

N  N  dz  N  dz 

Also,  the  following  sequences  of  functions  converge 
boundedly  to  the  indicated  limits  on  the  open  interval 
-H  <z  <0  as  N  +  ®: 

J.  -,2 

d  d  d  d 

—  gjz)  — g(z)  and  — :f  (z)  +  — rf(z)  .  (34) 

dz  “  dz  dz2  N  dz2 

(b)  Uniqueness  of  the  eigenfunction  expansions. 

A 

Given  a  sequence  of  complex  numbers  where 

n  =  ±1,  ±2,  and  so  on,  use  these  numbers  to  construct 

A  /V 

two  sequences  of  functions  {f..(z)}  and  (gv(z)}  by  means 

N  N 

of  the  symmetric  partial  sums 


(31) 


(32) 


(35) 


N 


fN(z)  =  l  Vn(Z) 
n  =  -N 


N 


gN(z)  =  l  Snikn*n(2) 
n  =  -N 


(36) 


Suppose  that  as  N  +  «>, 


f  (-H)  ->  f(-H)  , 

N 

r°  ^ 

fj^z)  h ( z)  dz  -+ 

-H 

r° 

g  (z)  h(z)  dz  -*• 
N 

J-H 


,0 

f  { z)  h (z)  dz  , 

-H 

f° 

g  (z)  h (z)  dz  . 
-H 


(37) 

(38) 


(39) 


whenever  the  test  function  h(z)  has  a  continuous  derivative. 

Then  a  =  for  all  n.  Indeed,  if  a  =  a  for  all  n,  then 
n  n  n  n 

conditions  (37-39)  are  satisfied  since  fN^z)  f  (z)  and 

gN(z)  -*■  g(z)  uniformly  on  the  closed  interval.  These 
convergence  conditions  therefore  completely  characterize 
the  eigenfunction  expansions. 

The  statements  in  the  existence  part  of  the  theorem 
will  be  proved  in  the  next  report.  The  uniqueness  of  the 
eigenfunction  expansions  is  easy  to  prove,  however,  and 
the  proof  motivates  formula  (30)  for  the  expansion 
coefficients.  Thus,  assuming  that  conditions  (37-39) 
are  satisfied,  we  base  this  proof  on  consideration  of 
the  following  quantity: 
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fO 

[k  f  (z)  -  i  g  (z)]  <j>  (z)P(~-}-dz  -  Hif  (-H)<fr  J-H)  .  (40) 
m  N  N  m  p  N  m 


It  is  clear  from  definitions  (35-36)  that  this  equals  the 
partial  sum 


N 

l  i 

n  =  -N 


n 


(k  +  k) 
m  n 


P(-H) 

<p(z)  4>  (z) - dz 


rm  n 


-H 


P  (z) 


-  H  i  <p  (-H)  4>_( -H) 
m  “ 


(41) 


In  view  of  relation  (19)  ,  only  the  n  =  m  term  can  contrib¬ 
ute  anything  nonzero  to  this  partial  sum.  Thus,  when  m  is 
fixed  the  limit  of  the  partial  sum  as  N  4  ”  is  equal  to  the 
value  of  this  term.  Also,  using  condition  (37)  and  taking 
h(z)  =  4>m(z)/p(z)  in  conditions  (38-39),  we  find  that  the 
limit  of  quantity  (40)  as  N  -*•  °°  can  be  expressed  in  the 
same  form  as  this  quantity,  but  with  the  functions  fN  and 
gN  replaced  by  f  and  g,  respectively.  Equating  these 
limits,  we  find  that 


-H 


p  (-H) 

(k  f  (z)  -  i  g(z)j  4>  (z) - dz  -  H  i  f  (-H)  4  (-H) 

m  m  p(z)  m 


=  a_ 


m 


f  0 


2  k 


m 


o  P  ( —  H)  .  9 

<fC(Z) - dz  -  H  i  4> “£( — H ) 

m  _ ,  m 


J-H 


P  (z) 


(42) 


Comparison  of  this  result  with  formula  (30)  shows  that 
am  =  am  as  claimed. 
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3.  FORMULATION  OF  THE  WAVEGUIDE  PROBLEM 

We  now  begin  to  consider  the  acoustics  problem  out¬ 
lined  in  the  introduction.  The  objective  in  this  problem 
is  to  find  how  a  waveguide  filled  with  fluid  responds  to 
small  amplitude  excitations  from  a  source  in  its  interior. 

We  let  the  top  boundary  of  the  waveguide  be  a  free 
surface  that  coincides  in  equilibrium  with  a  horizontal 
plane  of  infinite  extent.  One  direction  in  this  plane  is 
called  range  and  denoted  by  the  parameter  x.  All  features 
of  the  problem  are  translationally  invariant  along  the 
horizontal  direction  perpendicular  to  the  range.  The 
vertical  direction  is  denoted  by  the  parameter  z,  and 
z  =  0  on  the  top  boundary. 

The  bottom  boundary  of  the  waveguide  is  rigid  and 
uneven.  The  depth  of  the  bottom  varies  with  range  accord¬ 
ing  to  an  equation  of  the  form  z  =  -H(x),  where  the  depth 
function  H(x)  is  strictly  positive  and  has  a  continuous 
second  derivative. 

The  fluid  in  the  waveguide  is  inhomogeneous  but 
lossless.  The  local  mass  density  of  the  fluid  varies 
only  with  depth,  and  it  is  described  by  the  function  p(z). 
The  local  sound  speed  in  the  fluid  may  vary  with  range  as 
well  as  with  depth,  however,  and  it  is  described  by  the 
function  c(x,z).  The  second  derivative  of  p(z)  and  the 
first  partial  derivatives  of  c(x,z)  are  continuous. 

The  source  in  this  problem  is  concentrated  along  a 
line  that  is  parallel  to  the  direction  of  translational 
invariance.  It  generates  excitations  that  are  harmonic 
in  time.  Therefore,  the  pressure  response  of  the  wave¬ 
guide  can  be  expressed  as  the  real  part  of  the  complex 
product  exp(-iu>t)p(x,z)  ,  where  the  positive  constant  w 
is  the  angular  frequency  of  the  source  excitation,  and  the 
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parameter  t  is  time.  A  well-known  partial  differential 
equation  relates  the  response  function  p(x,z)  to  the 
excitation.  For  all  x  and  for  all  z  between  -H(x)  and  0, 
this  equation  has  the  form 


P(x,z) 


+  p (z) — 
3z 


1  9  . 

-  —  P (x, z) 

p (z)  3z 


+  K  (x,z)  p (x i z ) 


=  pg  6  (x  -  xg)  6  (z  -  zg)  , 


(43) 


where  K(x,z)  =  o)/c(x,z).  The  constant  p  indicates  the 
strength  of  the  source,  and  the  coordinates  x  =  xg  and 
z  =  zg  specify  its  position. 

The  pressure  response  is  subject  to  three  kinds  of 
boundary  conditions  in  this  problem.  One  of  the  conditions 
derives  from  the  vanishing  of  the  normal  component  of  fluid 
velocity  at  the  rigid  bottom.  Because  the  fluid  velocity 
is  proportional  to  the  pressure  gradient,  the  normal 
derivative  of  the  pressure  response  must  vanish  at  this 
boundary.  Vectors  that  are  normal  to  the  bottom  at  the 
position  z  =  -H(x)  are  parallel  to  the  vector  (H(x),l), 
where  H(x)  =  dH(x)/dx,  so  this  condition  can  be  formulated 
as 


3  •  3 

— p (x,-H (x) )  +  H ( x) — p (x , -H (x) )  =  0  (44) 

3z  3x 

for  all  x.  The  pressure  response  vanishes  at  the  top 
boundary,  and  we  note  that  its  tangential  derivative 
vanishes  there  as  well. 

p(x, 0)  =  0  ,  (45) 
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3 

—  p  (x,  0)  =0  , 
3x 


(46) 


for  all  x.  Finally,  the  pressure  response  far  from  the 
source  must  consist  essentially  of  outgoing  waves.  A 
mathematical  formulation  of  this  condition  appears  in  the 
last  section  of  this  report.  It  requires  that  the  functions 
H (x)  and  c(x,z)  converge  smoothly  to  finite  limits  as  |x| 
tends  to  infinity. 

H  (x)  -*•  H  ( ±°°)  and  H  (x)  -*•  0  as  x  -*■  ±°°, 

3 

c(x,z)  ■+  c(±°°,z)  and  — c(x,z)  ■+•  0 

3x 

uniformly  in  z  as  x  -*■  ±<*>  . 

4.  ASYMPTOTIC  LOCAL  ANALYSIS 

The  eigenvalue  problem  introduced  earlier  originates 
in  an  asymptotic  analysis  of  self-sustaining  oscillations 
in  the  waveguide.  This  analysis  produces  the  functions 
that  we  use  as  a  basis  for  expanding  the  exact  pressure 
response,  and  it  also  suggests  the  form  that  the  expansion 
takes.  It  is  of  interest,  therefore,  to  give  a  brief 
outline  of  the  asymptotic  analysis  and  to  discuss  some 
aspects  of  the  expansion  functions  in  this  context. 

The  aim  of  the  analysis  is  to  separate  the  variations 
in  the  pressure  response  that  depend  on  the  range  and  on 
the  depth.  In  most  cases  this  can  be  done  only  partially, 
so  we  adopt  the  following  representation  of  the  pressure 
response. 

p(x,z)  =  a  (x)  <Mx,z)  .  (49) 


(47) 


(48) 
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The  factor  a(x)  is  intended  to  isolate  the  major  range 
dependent  variations  in  p(x,z).  We  define  the  local 
wavenumber  k(x)  in  terms  of  this  factor. 


a  (x)  =  i  k (x)  a (x)  ,  (50) 

where  a(x)  =  da(x)/dx.  We  assume  that  k(x)  and  4>(x,z) 
are  slowly  varying  functions  of  range,  and  in  the  first 
stage  of  this  analysis  we  ignore  derivatives  of  these 
functions  with  respect  to  x  when  differentiating  p(x,z). 

For  example, 


3 

—  p  ( x ,  z )  =  a  (x)  i  k  (x)  4>(x,z)  , 

3x 


d  2 

— r;p(x,z)  =  -a(x)  k  (x)  4>(x,z)  . 

3x^ 


(51) 


(52) 


These  assumptions  turn  out  to  be  justifiable  whenever 
|k(x) |  H(x)  »  1  and  |H(x) [  «  1,  although  they  may  be 
valid  under  weaker  conditions  as  well. 

In  terms  of  the  functions  k(x)  and  <j>(x,z),  eqs.  (43-45) 
take  the  following  approximate  form  for  self-sustaining 
oscillations  in  the  waveguide,  that  is,  for  pg  =  0. 


,  ,  3 

' 

1 

3 

2  2 

P  (z)  — 
3z 

.  P  (z) 

—  <J>  (x,  z) 

3z 

+ 

K^(x,z)  -  K  (x) 

. 

=  0  , 

—  0 <x,-H(x) )  +  H (x)  i  k(x)  $ (x, -H (x) )  =  0  , 
3  z 

<Mx,0)  =  0  . 


(53) 

(54) 

(55) 
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Because  eqs. (53-55)  have  the  same  form  as  problem  (1-3) , 
the  set  of  solutions  k(x)  =  kn(x)  and  <Mx,z)  =  <pn(x,z)  is 
denumerable,  and  for  each  x  the  eigenfunctions  4>n(x,z) 
are  complete  in  the  two-fold  sense  described  in  the  prin¬ 
cipal  expansion  theorem.  Equations  (53-55)  completely 
determine  the  functions  kn(x) .  These  functions  are 
continuously  differentiable  and  slowly  varying  in  that 

kn(x)  'v  -  kn(x)  H  (x)  /  H  (x)  as  |n|  -*■  »  ,  (56) 

where  kn(x)  =  dkn(x)/dx. 

Problem  (53-55)  does  not  determine  the  dependence 
of  the  eigenfunctions  <J>n(x,z)  on  range,  however.  The 
next  stage  of  the  asymptotic  analysis  provides  this 
information.  We  set  k(x)  =  kn(x)  in  eq.  (50)  and  then 
attempt  to  find  a  small  and  slowly  varying  correction 
An(x,z)  that  will  reduce  the  difference  between  the  exact 
factor  <p  ( x ,  z )  in  eq.  (49)  and  the  approximation  4>n(x,z). 

<Mx,z)  =  $n(x,z)  +  An(x,z)  ,  (57) 

I  An(x,z)  |  «  |  <pn(x,z)  |  .  (58) 

When  we  differentiate  p(x,z)  at  this  stage,  we  retain  the 
first  derivatives  of  kn(x)  and  <f>n(x,z)  with  respect  to  x, 
but  we  ignore  the  derivatives  of  An(x,z)  with  respect  to 
x.  This  results  in  the  following  revisions  to  approxi¬ 
mations  ( 51-52 ) . 


—  p(x,z) 
9x 


a  (x) 


j^i  kn(x)  [$n(x,z)  +An(x,z)] 


_3_ 

9x 


4>n(x,z) 


9 
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+ 


(59) 


T 


b 


I 


ax" 


p(x,z)  =  a  (x) 


-  k^(x)  [<t>n(x,z)  +  An(x,z)] 


+  2  i  kjx)  —  <J>  (x,z)  +  i  k  (x)  <J>  (x,z) 
n  3x  n  n  n 


•  (60) 


By  substituting  these  approximations  in  eqs. (43-45) ,  we 
obtain  a  nonhomogeneous  boundary  value  problem  for  each 
An(x,z)  . 


a 

2,  .  .  2  . 

P  (z)  — 

A  (x, z) 

+ 

K  (x,z)  -  k  (x) 

3z 

p (z)  3z 

n 

*  „ 

=  -  2  i  k  (x)  —  A  (x,z)  -  i  k  (x)  <J> (x,z)  ,  (61) 

ax  n  n 

—  A  (x,-H  (x)  )  +  H  (x)  i  k  (x)  A  (x ,  -H  (x)  ) 

3z  n  n  n 

=  -  H(x)  —  <p  (x,-H(x)  )  ,  (62) 

3x  n 

A  (x , 0 )  =  0  .  (63) 

n 

In  the  present  discussion  we  are  interested  in  the 
condition  that  determines  whether  eqs. (61-63)  have 
solutions,  but  not  in  the  solutions  themselves.  This 
condition  is  a  corollary  of  the  theorem  in  appendix  A; 
it  states  that  a  particular  solution  of  eqs. (61-63) 
exists  if  and  only  if 


r  0 


_a_ 

3x 


( 


k  (x)  <f„(x,z) 
n  n 


p ( -H ( x) ) 
P  (z) 


dz  -  H (x)  —  —  4>2(x,-H(x)  ) 

2  ax  n 


-H(x) 
=  0  . 


J 
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(64) 


T 


Equation  (64)  is  equivalent  to  a  first  order  ordinary 
differential  equation  for  the  value  <t>n(x,-H  (x) )  on  the 
bottom  boundary,  and  it  determines  the  range  dependence 
of  each  eigenfunction  through  this  quantity.  This 
condition  also  causes  the  eigenfunctions  to  be  slowly 
varying  functions  of  range  in  the  sense  that 


—  <p (x,z)  =  max  |  <j>  (x,z)  |  ofk  (x)  H(x) 
3x  n  -H(x)<  z  <  0  l  n 


(65) 


Under  the  assumptions  of  the  asymptotic  analysis, 

the  oscillations  described  by  the  functions  k  (x)  and 

n 

<(>n(x,z)  are  independent  for  distinct  indices  n.  This 
property  motivates  us  to  use  these  functions  as  the 
basis  of  a  representation  of  the  exact  pressure 
response  in  the  waveguide  problem.  We  expect  the 
components  of  such  a  representation  to  interact 
weakly  on  the  whole. 

5.  EXACT  GLOBAL  ANALYSIS 

In  this  section  we  use  the  collection  of  solutions 
{kn(x)}  and  {4>n(x,z)}  of  problem  (53-55)  to  construct  an 
exact  representation  of  the  pressure  response  in  the 
waveguide  problem.  The  form  of  the  representation  is 
suggested  by  approximations  (49)  and  (51)  to  p(x,z)  and 
3p(x,z)/3x  in  the  preceding  analysis.  The  mathematical 
basis  for  this  representation  is  the  principal  eigen¬ 
function  expansion  theorem.  The  coefficients  in  the 
representation  turn  out  to  satisfy  a  linear  system  of 
first  order  ordinary  differential  equations,  and  they 
have  jump  discontinuities  at  the  position  of  the  line 
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source.  We  can  integrate  these  equations  to  find  the 
coefficients  from  initial  values  that  describe  the 
asymptotic  behavior  of  p(x,z)  far  from  the  source. 

In  order  to  develop  the  expansions  and  the  properties 
of  the  coefficients,  we  first  regard  the  pressure  response 
p(x,z)  as  known.  We  let  the  function  q(x,z)  denote  its 
partial  derivative  with  respect  to  range. 


—  p  ( x ,  z )  =  q  (x,  z)  .  (6( 

3x 

This  definition  makes  eq. (43)  formally  first  order  in  the 
operator  3/3x. 


—  q(x,z)  +  p(z)  — - - - — —  p  ( x ,  z )  +  K2(x,z)  p  (x  ,  z) 

3x  3z [p (z)  3z 

=  p  6 (x  -  x  )  6 (z  -  z  )  .  (67) 

s  s  s 

For  fixed  x,  the  boundary  conditions  (44-46)  imply 

that  the  functions  f(z)  =  p(x,z)  and  g(z)  =  q(x,z)  satisfy 

requirements  (27-29)  of  the  principal  eigenfunction 

3  3 

expansion  theorem.  The  derivatives  3  p(x,z)/3z  and 
2  2 

3  q(x,z)/3z  are  continuous,  except  at  the  source,  because 

2  2 

the  derivatives  3c(x,z)/3x,  3c(x,z)/3z  and  d  p(z)/dz  are 
continuous.  Therefore,  formula  (30)  uniquely  determines 
the  coefficients  an(x)  in  the  following  expansions. 


p(x,z)  =  l  a^x)  <J>n(x,z)  , 

n  =  -oo 


q  ( x ,  z )  =  l  an(x)  i  k  (x)  4>n(x,z)  . 
n  =  -oo 
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The  expansion  of  the  depth  component  of  the  pressure 
gradient  is  obtained  from  expansion  (68)  by  termwise 
differentiation. 


—  p(x,z)  =  l  a  (x)  —  <j>  (x ,  z )  .  (70) 

n  =  -oo  n  ^  n 

Expansions  (68-70)  stand  for  the  limits  as  N  ■*  ®  of  the 
corresponding  symmetric  partial  sums  from  n  = -N  to  n=N, 
and  in  this  sense  they  converge  uniformly  on  the  closed 
interval  -H(x)  £  z  <_  0  if  x  ^  xg.  Expansions  (69-70)  may 
be  regarded  almost  as  representations  of  the  velocity 
field  in  the  fluid  since  the  fluid  velocity  is  propor¬ 
tional  to  the  pressure  gradient. 

We  can  apply  formula  (30)  to  find  an  expression  for 

each  difference  a  (x)  -  ik  (x)  a  (x)  ,  where  a  (x)  =  da  (x)/dx. 

n  n  n  n  n 

Let  us  define  two  auxiliary  functions  f(x,z)  and  g(x,z) 
in  terms  of  these  differences. 


f(x,z)  =  l  (an(x)  -  i  kn(x)  an(x)]  4>n(x,  z)  ,  (71) 

n  = 

00 

g(x,z)  =  l  [an(x)  -  i  kn(x)  an(x)]  i  kn(x)  4>n(x,  z)  .  (72) 

n  =  -oo 

These  functions  are  constructed  like  expansions  (25-26) , 
and  their  convergence  properties  are  in  accord  with 
conditions  (37-39)  in  the  uniqueness  part  of  the  principal 
expansion  theorem.  This  is  due  to  the  asymptotic  relation 
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an(x)  ^  i  k_n( x )  an(x)  as  |n 


->  OO 


(73 


which  generalizes  eq.(50)  in  the  preceding  asymptotic 
analysis.  Consequently,  formula  (30)  implies  that 

an(x)  -  i  kn(x)  an(x)  =  Cn(x)/Dn(x)  ,  (74 


where 


cn<x)  = 


[k  (x)  f(x,z)  -  ig(x,z)]  <J)  (x,z)  —  -  dz 

n  n 


P  (z) 


-H(x) 


-  H  (x)  i  f  (x,-H  (x)  )  4>n(x,-H  (x)  )  , 


(75) 


Dn(x) 


2  kn(x) 


2  .  P(-H(x) ) 

t>n(x,z)  - dz 


P  (z) 


-H(x) 


-  H  (x)  i  <jT(x,-H(x)  )  . 


(76! 


We  use  eq. (74)  to  derive  the  linear  system  of  first 
order  ordinary  differential  equations  for  the  coefficients 
First,  we  construct  new  but  equivalent  expressions  for  the 
auxiliary  functions  f(x,z)  and  g(x,z)  by  substituting 
expansions  (68-69)  in  relations  (66-67)  and  rearranging 
the  resulting  sums.  The  indicated  substitutions  produce 
two  constraints  on  the  coefficients  and  their  derivatives. 
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l'  an(x)  4>n( x ,  z ) 
n  =  -<» 


+ 


n 


f  an(x) 


3 

—  4>n(  x ,  2 ) 
3x 


CO 

l  an(x)  i  kn(x)  4>n(x,  z)  , 
n  =  -<*> 


(77) 


oo  00 

l  a  (x)  i  kn(x)  4>n(x,z)  +  l  an(x)  i  yx)  4>n(x,z) 
n  =  -oo  n  =  -°° 


+  l  an(x)  i  kn(x)  — <j>n(x,z)  + 

3x 


n  =  -oo 


OO 

l  an(x)  k^(x)  4>n(x,z) 

=  — OO 


=  p  5  (x  -  x  )  6  (z  -  z  )  . 
s  s  s 


(78) 


Therefore,  if  x  /  x  then  the  auxiliary  functions  defined 

s 

in  eqs. (71-72)  have  the  following  equivalent  representations. 


f(x,z)  =  -  l  a  (x)  — <p (x,z)  , 

„  _  „  3X  “ 

n  = 


(79) 


g(x,z)  =  - 


y  a  (x)  i 

L  n 


n  =  -°° 


kn(x)  *n(x,z) 


+  k  (x)  —  <()  (x, z) 
n  3x  n 


(80) 


Substituting  these  expressions  for  f(x,z)  and  g(x,z)  in 
eq.  (75)  yields  an  expansion  of  the  factor  Cn(x)  in  terms 
of  the  coefficients.  This  transforms  relation  (74)  into 
a  differential  equation  that  is  linear  in  the  coefficients. 
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If  x  f  xg  then 

00 

an(x)  -  ikn(x)an(x)  =  -  \  am(x)  cmn(x) /  Dn(x)  ,  (81) 

m  =  —oo 


where 


t  0 


c  (x) 
mn 


=  [k  (x)  +  k  (x)] 
1  m  n  1 


3  .  . 

—  <P  Jx,z) 
3x 

p  (-H(x)  ) 

<p  (x,z)  - dz 

p(z) 


-H(x) 


f  0 


+  km<x) 


p (-H (x) ) 

t>m(x,z)  <f>n(x,z)  - dz 


P  (z) 


-H(x) 


-  H (x)  i 


9  x 


<j>n(x,-H(x)  ) 


(82) 


for  m,n  =  +1,  ±2,  and  so  on. 

It  is  interesting  to  note  some  structural  features 
of  these  equations.  First,  eq. (81)  has  the  same  form 
as  standard  coupled  mode  equations  associated  with  wave¬ 
guides  that  have  parallel  plane  boundaries.  On  the  other 
hand,  eqs. (76)  and  (82)  differ  from  standard  expressions 
because  they  include  terms  that  are  proportional  to  H(x). 

In  the  derivation  of  these  equations  we  have  not  used 
condition  (64),  which  determines  the  range  dependence  of 
the  eigenfunctions.  If  this  condition  is  adopted,  however, 
it  causes  the  diagonal  matrix  elements  cnn(x)  to  vanish. 
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If  x  =  x  then  the  derivatives  an(x)  fail  to  exist 
because  each  coefficient  an(x)  has  a  jump  discontinuity 
at  this  range.  We  calculate  the  jumps  as  follows. 

an(xs+0)  “  an(xs-0) 

=  lim  [an(xg  +  e)  ~  an(xg  ~  e)l  •  (83) 
e  -*  0 

A  pair  of  constraints  on  these  jump  discontinuities  can 
be  obtained  by  integrating  eqs. (77-78)  from  x  =  xg  -  e  to 
x  =  xg  +  e  and  then  letting  e  0+  in  the  result. 


I  [a  (x  +  0)  -  a  (x  -  0)]  <j>  (x  ,z)  =  0  ,  (84) 

1  n  s  n  s  J  Yn  s 

n  =  -°° 

oo 

l'  (an(xs+0)  “  an(xs'0)]  ikn(xs}  ♦n(xs»z) 

n  =  -oo 

■  Pg  5 (z  "  zs)  .  (85) 

Because  these  constraints  have  the  same  form  as  eqs. (25-26) , 
we  can  evaluate  the  jump  discontinuities  by  another  applica¬ 
tion  of  formula  (30)  with  f(z)  =  0  and  g(z)  =  p  6(z-z  ). 

s  s 

an(xs+0)  '  an(xs'  0) 
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1<*>n(xs'zs)  P  (~H  (xs) ) 

-  -  Pg - .  (86) 

Dn(xs}  P(zs) 


Although  the  pressure  response  has  been  regarded  as 
known  up  to  this  point,  it  is  in  fact  initially  unknown 
in  the  waveguide  problem.  Reversing  the  direction  of  the 
preceding  discussion  leads  to  a  method  for  finding  p(x,z) 
in  this  problem.  First,  we  find  the  coefficients  an(x) 
by  integrating  eqs.(81)  subject  to  conditions  (86),  and 
then  we  use  the  coefficients  to  construct  expansion  (68) 
for  p(x,z).  In  order  to  start  the  integration,  however, 
we  need  values  of  the  coefficients  at  some  reference 
points.  Our  choices  of  the  initial  conditions  are 
motivated  by  the  requirement  that  the  pressure  response 
should  consist  essentially  of  outgoing  waves  far  from 
the  source. 

Under  assumptions  (47-48)  of  the  waveguide  problem, 
the  functions  an(x)  turn  out  to  be  basically  exponential 
far  from  the  source.  Each  kn(x)  tends  to  purely  real  or 
purely  imaginary  limits  kn(±<»)  as  x  -*  +“>,  and  there  are 
complex  amplitudes  An(±®)  such  that 

a  (x)  exp{-  ik„(±°°)  x}  •+  A  (±°°)  as  x  ■+■  ±°°  .  (87) 

n  *■  n  n 

Due  to  the  underlying  exp(-ioot)  time  dependence  of  the 
pressure  response,  the  following  initial  conditions  at 
the  indicated  ranges  exclude  propagating  waves  that 
approach  the  source  from  infinity.  These  conditions 
also  exclude  stationary  oscillations  that  grow  unbound¬ 
edly  in  amplitude  as  |x|  increases  without  limit. 

a  (-<*>)  =  0  if  Re{k  (-»)}  >  0  or  Im{k  (-»)}  >  0  , 
n  n  n 

a  (+<*>)  =  0  if  Re{k _(+<*>)}  <  0  or  Im{k  (+<*>)}  <  0  . 
n  n  n 


(88) 

(89) 
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APPENDIX  A:  A  NON HOMOGENEOUS  BOUNDARY  VALUE  PROBLEM 


In  this  appendix  we  state  a  necessary  and  sufficient 
condition  for  the  existence  of  a  particular  solution  <J>(z) 
to  the  following  nonhomogeneous  boundary  value  problem. 


d 

o(z)  _ 

1 

d 

—  <*>  ( z ) 

+ 

K2(z)  -  k2 

dz 

.  P  (z) 

dz 

n 

for  -H  <  z  <  0  , 

(Al) 

1$(-H)  +  Hik  4>(-H)  =  -  H  y  , 
dz  n 

(A2) 

<t  ( 0 )  =  0  . 

(A3) 

In  this  problem,  kn  is  an  eigenvalue  of  problem  (1-3). 

We  assume  that  the  function  F(z)  has  a  continuous  second 
derivative  on  the  closed  interval  -H  <_  z  £  0  and  that  F(0) 
vanishes.  The  parameter  y  may  take  any  complex  value. 

Theorem. 

Problem  (A1-A3)  has  a  solution  if  and  only  if 

,0 

F(z)<Uz)  dz  -  Hy4(-H)  =  0  ,  (A4) 

P  (2) 


where  the  eigenfunction  4>n(z)  is  associated  with  kn>  In 
addition,  if  relation  (20)  holds  for  all  indices,  then 
one  particular  solution  can  be  represented  as  follows. 


oo 

♦<*>  -  r  •  <“> 

m  =  -*> 

(m  /  n) 
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I 


p(-H) 


1 0 


F  ( z )  4>m(  2  ) 


-H 


P  (z) 


dz  - 


H  y  <J>  ( -H) 
m 


am 


2  k 


m 


2  P(-H) 

<T(z)  -  dz 

m 


-H 


P  (z) 


H  i  02(-H) 
m 


(A6) 


It  is  easy  to  prove  the  necessity  of  condition  (A4)  . 
As  in  the  derivation  of  relation  (19) ,  we  first  express 
eq.  (1)  for  *f>  ( z)  and  eq.  (Al)  for  $(z)  in  the  respective 
operator  forms  M4>n,kn)  =  0  and  L(4>,kn)  -  F  =  0.  Then 
we  obtain  the  following  identity  by  partial  integration 
of  the  difference  ($/P )  L  ( 4>n,kn)  -  (<j>n/p)  L(<J>,kn)  : 


f° 

F(z)  <f»  (z)  -i—  dz 
n  . 


1  d 

<t>(  z ) - <t  ( z ) 

n  p  ( z )  dz 


1  d 

*(z) - <tn(z) 

p  (z)  dz 


z  =  0 


z  =  -H 


(A7) 


This  identity  leads  directly  to  condition  (A4)  after  we 

apply  boundary  conditions  (2-3)  for  <{>n(z)  and  boundary 

conditions  (A2-A3)  for  <t(z). 

One  can  prove  the  sufficiency  of  condition  (A4)  by 

verifying  that  eqs. (A5-A6)  represent  a  particular  solution 

of  problem  (A1-A3)  when  it  holds.  This  proof  requires  use 

of  the  principal  eigenfunction  expansion  theorem,  but  it 

will  not  be  given  here.  Equations  (A5-A6)  are  provided 

only  to  indicate  the  form  of  the  correction  terms  A  (x,z) 

n 

in  the  asymptotic  analysis  of  the  waveguide  problem. 
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